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Zagadnienie optymalizacji jonowymiennego rozdziału pierwiastków ziem 
rzadkich 


Проблема оптимализации ионообменного разделения редкоземельных 
элементов 


The Problem of Optimization of Ion Exchange Separation of Rare Earth Elements 


Na jakość rozdziału jonowymiennego wpływają dwie grupy czynni- 
ków: a) czynniki określające zdolności rozdzielcze mieszaniny (charakte- 
ryzuje je współczynnik separacji); b) czynniki określające sprawność 
kolumny, którą można charakteryzować za pomocą EHTP (efektywnej 
wysokości teoretycznej półki). Sprawność kolumny jest tym lepsza, im 
mniejszą wartość EHTP ma dana kolumna. 

Pomysł półki teoretycznej został zaczerpnięty z teorii destylacji 
i wprowadzony po raz pierwszy do chromatografii przez Martina 
i Syngea [1] w r. 1941 (teoria ich nie jest obecnie używana przy 
dyskusji rozdziałów jonowymiennych i z tego powodu nie będzie tu 
omawiana). Dalszych modyfikacji dokonali Mayer i Tompkins 
[2, 3], stosując teorię półek do dyskusji rozdziału pierwiastków ziem 
rzadkich na kolumnach jonowymiennych. Oparli się oni na uproszczo- 
nym założeniu, przyjmując nieciągłość przepływu, przez co, szczególnie 
przy niskich wartościach współczynnika podziału, teoria ich nie daje 
dobrej zgodności pomiędzy parametrami wyliczonymi i otrzymanymi 
na drodze eksperymentalnej [4]. Glueckauf [4, 5] rozszerzył 
i uogólnił teorię półek, wprowadzając zmienne ciągłe, przez co model 
stał się bardziej adekwatny do rzeczywistości. 

Teorię półek przedstawiono w następującym porządku: 

1. Określono wysokość EHTP na drodze teoretycznej w oparciu 
o podstawowe wielkości, tj. współczynniki dyfuzji, rozmiary ziarna jo- 
nitu, szybkość przepływu eluentu i współczynnik podziału. Otrzymane 
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zależności umożliwiają badanie wpływu wymienionych wyżej parametrów 
na przebieg procesu. 

2. Podano teorię krzywych przebicia, która pozwala powiązać EHTP 
z rozkładem rozdzielanej substancji na kolumnie, a także daje możliwość 


obliczenia 


wysokości półki z danych eksperymentalnych. W oparciu o wy- 


prowadzone relacje można badać wpływ współczynnika separacji, współ- 


czynnika 
działu. 


podziału i EHTP na kształt krzywej przebicia i czystość roz- 


3. Określono długość kolumny dla przeprowadzenia rozdziału z uzys- 
kaniem produktów o założonej z góry czystości przy danym współczyn- 
niku separacji. 


Symbolika: 

$ — przekrój poprzeczny kolumny (cm?); 
l  — liniowa długość kolumny (ст); 

X — objętościowa długość kolumny (ст); 


Vie = 


У 
p = 


długość pasma zaabsorbowanej substancji na wierzchołku 
kolumny; 
rozmiary ziarna jonitu; 
liczba półek teoretycznych, odpowiadająca danej kolumnie; 
— efektywna wysokość teoretycznej półki; 
objętość kolumny odpowiadająca jednej półce teoretycznej 
АХ = Л: S); 
liczba półek teoretycznych dla eluentu; 
kolumnowy współczynnik podziału ( =% koráli 

і 
ilość substancji zawarta w 1 ст? kolumny (mval: ст- 3); 
stężenie substancji w roztworze (mval · cm”* roztworu); 
stężenie substancji w јопісіе (mval*cm"$ jonitu) — jest to 
ilość substancji zawarta w jonicie mieszczącym się w 1 ст? 
kolumny; 
stosunek roztworu zawartego między ziarnami jonitu w ko- 
lumnie do całkowitej objętości kolumny; 
objętość jaka potrzebna jest do wymycia punktu o stężeniu 
с =: 
pojemność jonitu (туа!|* ст-3 jonitu); 
ułamek molowy substancji і zawartej w roztworze; 
współczynnik podziału Mayera i Tompkinsa; 
współczynnik separacji; 
szerokość pasma elucyjnego w punkcie С; = 0,368 Cmax: (cm3); 
objętość wymycia piku (ст); 
współczynniki dyfuzji: w roztworze D; oraz w fazie jonitu Он; 
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5 — grubość warstewki adhezyjnej; 

Е — objętościowa szybkość przepływu eluentu (ст? · зек-1; 

Е — liniowa szybkość przepływu eluentu (cm · sek 1); 

С; — stężenie substancji w roztworze w równowadze z jonitem 
(mval + ст); 5 

а; — stężenie substancji w fazie jonitu w równowadze z roztworem 
(mval · cm~? jonitu); 

V; — objętość wycieku (стз); 

гп — ilość substancji zawarta w fazie jonitu na półce n przy porcji 


roztworu T; 


Ln — ilość substancji zawarta w roztworze na półce п przy porcji 
roztworu T; 


a — ułamek charakteryzujący obciążenie kolumny (= Э; 
: rób z : aż N* 
р — parametr charakteryzujący obciążenie kolumny = 2 |: 


— szybkość migracji pasma wzdłuż kolumny. 


RÓWNANIE BILANSU DLA ZMIENNYCH CIĄGŁYCH [4, 5, 6] 


W dowolnym elemencie kolumny Ax (ryc. 1) zawarta jest ilość sub- 
stancji określona przez: 


ДАСТОР =т= Ту 


Ryc. 1. Schemat kolumny 


Przepłynięcie ilości roztworu dv spowoduje, że przez powierzchnię 
przekroju kolumny w punkcie х-Ах wpłynie do tej warstewki ilość 
substancji określona przez: 
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о о тылзаш kT бы (2) 

Ilość zaś substancji, jaka wypłynie w tym samym czasie z warstewki Ax 
przez powierzchnię przekroju kolumny w punkcie x, wyniesie: 

ФС зле шалы аай солод) 

Przenoszenie substancji wzdłuż kolumny zachodzi ponadto na drodze 

dyfuzji. W czasie dt, jaki jest potrzebny na to, aby ilość roztworu dv 


przeszła przez warstewkę Ax, zostanie dostarczona na drodze dyfuzji 
do tej warstewki ilość substancji określona przez: 


D,:s:f:dt | 2Сх—Ах) 
ma w F idia (4) 
W tym czasie wypłynie z niej ilość określona przez: 
D,*s:f:dt [3C (2) А 
s (Be), ш) 


Przepływy te spowodują zmianę zawartości substancji w Ax, której 
ilość po przejściu dv roztworu wyniesie: 
Кс а т % |... (6) 


Z powyższych rozważań wynika, że zmiana (5Q;) zawartości substancji 
w pewnym segmencie kolumny Ax jest spowodowana przez przepływa- 
jący roztwór (OQ;)przept. oraz przez dyfuzję wzdłuż kolumny (SQ;)ayr., 
co można zapisać: 


(56) = МӘ а. BPO гаса. ан ае. (® 
Uwzględniając równania (1—6), znajdujemy: 
(29) = Ах [Q, (v+dv) — ©, (v)] (7a) 
(EQ;)przepr. =dv [С, (ж—Ах) a g С, (=)] (7b) 
о), S D,'s-f:dt | | н A | OC, 4 
y2 дх » дк |, (7с) 


Po rozwinięciu wyrażeń (7а, с) w szereg Taylora, wstawieniu do 
równania (7) i przeniesieniu wszystkich wyrazów na jedną stronę, otrzy- 
mamy: 


ь у (2509 НО 
n! дъ z n! дх v 


D; f-dv >| ә” СЫ) | 
ЕА Е —А FI — 0 
Ру? 9х" ( z) 
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gdzie: wre s= 


Ostatnie. równanie jest i niewygodne w zastosowaniach praktycznych. 
Po odrzuceniu w nim dalszych składników szeregu Taylora, jako war- 
tości bardzo małych, przyjmie ono“znacznie prostszą postać [5, 6, 7]: 


© ©, (0) ш ӘС,(а)\ _ Ах / 9°С,(®)\ _ DfAx | O*C,(x)| _ 
ойлаша чыте | mal 28 |7 


#5 ‚.(8) 


2 


OKREŚLANIE ЕНТР NA DRODZE TEORETYCZNEJ 


Równanie EHTP dla elucji. W trakcie przepływu roz- 
tworu przez kolumnę jego czas kontaktu z ziarnem jonitu jest ograni- 
czony i z reguły zbyt mały na to, aby zdołała ustalić się równowaga. 
Odchylenia od stanu równowagi będą zdeterminowane przez najpowol- 
niejszy z procesów warunkujących wymianę jonową. Będą to: dyfuzja 
wewnątrzcząsteczkowa lub dyfuzja w warstewce adhezyjnej, albo jedna 
i druga jednocześnie [6—15]. 

Najpierw zajmiemy się przypadkiem, gdy główną rolę odgrywa dy- 
fuzja wewnątrzcząsteczkowa. W równaniu bilansu [8] występuje wyra- 


żenie 2 9. fe» | „natomiast szybkość dyfuzji zwykle wyraża się przez 
z 


ИЛ! 


Dlatego też przy wprowadzaniu równania dyfuzji do równania bilansu 
konieczne jest znalezienie relacji typu: 


Сө» ҮГӘ, o 


Jak wykazał Glueckauf [5, 15], z dobrym przybliżeniem można 
stosować dla dyfuzji wewnątrzcząsteczkowej następujące równanie: 


- © (0) Dri |. 
| ЕГ | oz 40-00] 


Równanie to po uwzględnieniu zależności: 


а, (®) = Qy(v) — ЈС (х) 
qio) =K4C;(2) 


(9b) 


daje się przekształcić w relację typu (9), co po wstawieniu do równania 
bilansu (8) daje: 
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ӘС (а) zeig) 0,071 Fr? 
| 3 } + (Кс | m a BARK +fF Е D,, | 4 
элн үн, УУ 


2?C,(2) _ Ах 3°C (x) kd рјАх ӘС) = (10) 
к | да? |, zł dx? |. w | да? | 


Natomiast dla dyfuzji w warstewce adhezyjnej można zastosować [5] 
równanie: 


да 
= c 
| д е} p — a _ (9е) 
Równanie to ро uwzględnieniu relacji: 
NBL TA КАР ш” 
1707р" 


qv) = ©(0)+ fC,(z) 


pozwala па określenie zależności typu (9). Wstawiając następnie do rów- 
nania bilansu (8) otrzymamy: 


д Сух) ə C,(v) Кс’ 0,133r,*F 3C (x) 
En Ох |, ЫК Е дъ |, шз sam azja ox” р 


звы 
i == e (55) Żar an 
2 | да? |), Fy2 да? |, 


Równania (10) i (11) różnią się czynnikami zaznaczonymi klamrami, co 
sugeruje, że czynniki te określają wpływ na stan równowagi odpowiednio 
dyfuzji wewnątrzcząsteczkowej i dyfuzji w warstewce adhezyjnej. Stąd 
dla procesu, w którym odgrywają one równorzędną rolę, równanie bi- 


lansu przyjmie postać: 
СД 


EHTP, 
———L-—— 
SW). эзше эрг) s | э) Кс 0.142r*F 
| 9 (x) 5A ar dv |, 2\ дз |). ‘| EE E | 
ЕНТР, EHTP,  EHTP, 
p ZZ ——— > = 
‚8: 0,266r,* Е V2 Di f 
— A: M ——Ę+ M, + jo (12) 
Е (Kot f) тЫ „+ Е 


Wyrażenia wzięte w nawias kwadratowy określają wpływ na EHTP 
poszczególnych czynników, tj. dyfuzji wewnątrzcząsteczkowej (ЕНТРЬ), 
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dyfuzji. w warstewce adhezyjnej (EHTP;), rozmiarów ziarna jonitu 
(ЕНТР.), dyfuzji wzdłuż kolumny (ЕНТР,), gdzie: 


Ка 0,142т Е 


EHTP о ———— 3 
"ЕТ D, ыы 
Кс’ 0,2667 F 
ЕНТР, =——— l (14) 
(Ko + f)? D(1+70r F) 
ЕНТР == 1,64т, (15) 
Dj:V2 
EHTP, = Е (16) 


Równanie EHTP dla rugowania. Rugowanie następuje 
wówczas, gdy jon wypierający ma większe powinowactwo do jonitu od 
jonu wypieranego. Wtedy po pewnym czasie, zależnym od współczyn- 
nika separacji, ustala się stacjonarna izochrona [6, 8, 16—24]. W stanie 
stacjonarnym w wyodrębnionej warstewce Ax zawarta jest ilość substan- 
cji określona przez: 

Ах (qo FfC;) 


W jednostce czasu opuści Ах ilość substancji równa s-F' С, stąd czas 
potrzebny na jej przejście do następnej warstewki wyraża się przez: 


Аза, ус) _ 40а, +fC) 


At= = 
sFC, FC 


i dalej: 


к ООА, ЕС, | Ж Е 
=” AR ЕС ДЫЛ Ич 
а Р CEL 


W tym przypadku podział substancji pomiędzy fazę janitu i roztwór 
określony jest przez zależność qo/ci. Przy podobnym postępowaniu јак 
w równaniach (9—12) z zastąpieniem Ко przez qo/ci dochodzimy do 
określenia EHTP dla rugowania: 


w (17) 


q /C: 0,142r.?F 


EHTPp = (°С, +77 Dr, (18) 
а`/С | 0,266т„°Е 
ЧТ 0200, В 18) 
BEN күз +f) D;(1+70r,F) 
ЕНТР, = 1,64 r, (20) 
ЕНТР, = spya (21) 
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Uwagi ogólne. Jak wykazano wyżej (równania 13—16, 18—21) 
wartość ЕНТР (=EHTP, +ЕНТР;+ЕНТР,+ЕНТР,) zależy od: tempera- 
tury, szybkości przepływu eluentu, rozmiarów ziarna jonitu i współ- 
czynnika podziału. Pierwsze trzy czynniki wiążą się z szybkością wy- 
miany jonowej. Okazuje się, że wraz z przybliżaniem się reagującego 
układu do stanu równowagi EHTP maleje. 

Na kolumnie jonowymiennej istnieją zasadniczo dwie tendencje. 
Z jednej strony wszystkie procesy zdążają do stanu równowagi, z dru- 
giej — wraz z zachodzącym rozseparowywaniem substancji powstaje 
sukcesywnie wzdłuż kolumny gradient stężenia. Wywołany w taki sposób 
stan nierównowagowy będzie powodował wzdłuż kolumny dyfuzję 
(w roztworze zawartym pomiędzy ziarnami jonitu), której działanie jest 
niekorzystne dla rozdziału, gdyż prowadzi do rozmywania się izochrony. 
Na szczęście jest to proces na tyle powolny, że w normalnej tempera- 
turze można go, bez popełnienia większego błędu, pominąć. 

Wzrost wartości współczynnika dyfuzji wraz ze wzrostem tempera- 
tury prowadzi do zmniejszenia się EHTP dla jonów wielowartościowych 
(ryc. 2a i 6), gdyż obniżanie się wartości ЕНТР, і EHTP; (równania 13, 
14, 18, 19) zachodzi szybciej niż wzrost EHTP4 (równania 16, 21), jed- 
nakże przy ruchliwych jonach wzrost ЕНТР, przewyższa obniżanie war- 
tości ЕНТР, i EHTP; (ryc. 2b). 


12 
b a b 
10 
o 

эк. 
са 
26 0 6 
Е £ 
E 4 4 
= © > 
= к о 

2 о 2 102009077 

u 
о о Е 
о 20 «o бо во ' 100% о 20 40 60 во 100 + 


Ryc. 2. Efektywna wysokość teoretycznej półki jako funkcja temperatury; a — dla 
pierwiastków ziem rzadkich [27], b — dla potasowców [27] 


Zmniejszenie szybkości przepływu eluentu wydłuża jego czas kon- 
taktu z ziarnem jonitu, a tym samym przybliża reagujący układ do 
stanu równowagi, przez co powoduje obniżanie wartości EHTP (ryc. 3, 
4, 6), co oznacza polepszenie sprawności kolumny. Jednak obniżanie war- 
tości EHTP poprzez zmniejszanie szybkości przepływu eluentu prowadzi 
do wydłużenia czasu rozdziału, co z ekonomicznego punktu widzenia może 
okazać się nieopłacalne. 
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Ryc. 3. Efektywna wysokość teoretycznej półki jako funkcja szybkości przepływu 


eluentu (Е jest wyrażone w ml/min.) [26]; 1 — KG=0,6, 2 — Kc=30 


EHTP(wmm) 


ama 


o ол аз оз 
fo (w mm} 


Ryc. 4. Efektywna wysokość teoretycznej półki jako funkcja rozmiarów ziarna 


jonitu przy różnych szybkościach przepływu eluentu [25]; A — 5 ml/sek. cms, 
В — 2 ml/sek. cm? С — 1 ml/sek. ст?, D — 0,3 ml/sek. cm? Е — 0,1 ml/sek. cm? 
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Ryc. 5. Efektywna wysokość teoretycznej półki jako funkcja wartości współczynnika 
podziału; na osi rzędnych odłożono liczbę półek teoretycznych dla kolumny o dłu- 
gości 37 ст (N=37/EHTP) [26]; 1 — elucja Sr, 2 — elucja Sr, 3 — elucja Eu 


Zmniejszenie rozmiarów ziarna jonitu prowadzi do obniżenia wartości 
EHTP (ryc. 4). Jednak stosowanie jonitu o zbyt rozdrobnionym ziarnie 
stwarza duże trudności natury technicznej, a w krańcowym przypadku 
wręcz jest niemożliwe ze względu na wysoki opór, jaki stawia takie złoże 
przepływającemu roztworowi. 


Wzrost współczynnika podziału prowadzi do obniżenia wartości EHTP 
(туе. 5, 6). Polepszeniu sprawności kolumny w ten sposób towarzyszy 


= 827 
F 


wzrost czasu rozdziału (: tak że w końcu dopiero względy 


ekonomiczne konkretnego procesu zadecydują o optymalnym wyborze. 


Reasumując, można powiedzieć, że określenie EHTP w oparciu o pod- 
stawowe parametry zasadniczego procesu wymiany jonowej jest jedną 
z głównych zalet tej teorii, gdyż pozwala bezpośrednio badać wpływ 
tych parametrów na jakość rozdziału, a w połączeniu z rozważaniami 
natury ekonomicznej dobierać najbardziej opłacalne warunki do prowa- 
dzenia procesu. 
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EHTP х 10? ст 


о 20 4O 60 80 100 120 140KG 


EHTP х 10%cm 


о 20 40 60 во 100 120 10 Kg 

0-1 4-2 0-3 

v- x-5 0-6 
Ryc. 6. Efektywna wysokość teoretycznej półki jako funkcja współczynika podziału 
przy różnych wartościach temperatury i szybkości przepływu eluentu [28]; a — szyb- 
kość przepływu eluentu 1,8—2,1 :10—2 ml*sek.—1 cm—?; 1 — elucja У, La, Ce i Pr 
w temp. 18°С, 2 — elucja Y, La i Pr w temp. 50°С, 3 — elucja Tb, Gd i Eu 
w temp. 55°С, 4 — elucja Lu, Sc i Ho w temp. 65°С; b — temp. 25°С, 5 — Y, La 


i Ce przy F= 1,82 -10-2 ml: sek.—!1 cm—?2; 6 — elucja Tu, La i Tb przy F= 1,06 · 10-2 ml - 
ре Вы = 


TEORIA KRZYWYCH PRZEBICIA 


Teoria Glueckaufa dla elucji. Po odrzuceniu w rów- 
naniu bilansu (8) składnika odnoszącego się do dyfuzji wzdłuż kolumny 
oraz po uwzględnieniu: 


Q= (К. +) С, =4C, 


Хх, _ ы У, тох дү 
Anadia госте $ 
równanie bilansu przyjmie postać [29]: 
oC, С OC 
i + OC, b, a Ay doj (22) 
ЭМ, ] №, дм, м, 2 \9М, ]м, 


Załóżmy, że kolumna zawiera na swoim wierzchołku pasmo zaabsor- 
bowanej substancji o długości №. Wtedy warunki początkowe i i graniczne 
będą miały postać [24]: 
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C(ON, <N)=C; 

c,(O,N, > N)=0 

c,(M„O) =0 
Rozwiązując równanie (22) przy wyżej wymienionych warunkach, otrzy- 
mamy: 


с Ni— Mi Ni Ni — Mi 
T „A VM Ае Т _ 


gdzie: 
A = Г + d 
= 2 mz» (z) t 


Wartości funkcji A {t} są stabelaryzowane [30, 31, 32]. Dla N, ŚNi 
równanie (23) można doprowadzić [29] do prostszej postaci: 


пи /N; N: (Wi—V3)? 
а= |у -3 = 24 
i V: ] 2% р | 2 уу, (24) 
Równanie to osiąga maksimum, gdy exp (0)=1 i wtedy: 
Le тч м 
С: = Cmax у. y= 
stąd: 
ту \2 
Е. ре (25) 
m 


Podstawiając do równania (24) C;= 0,368 Cmax, a pęsiępnię logarytmując 
otrzymamy: 


__ ЗУ У 8 
(Vi— Vie)? № 

Uproszczenia natury matematycznej, wprowadzone przy przejściu od 

równania (23) do równania (26), u których podstaw leżało założenie, że 

№,°<< Ni ogranicza zastosowanie równania (25) i (26) do przypadków, gdy 

stosunek części obciążonej kolumny do części rozwijającej jest mały. 
Wielkość jego jest ograniczona przez [29]: 
M лы 

ИЗМ! 

W późniejszych badaniach Glueckauf [33] rozszerzył zastosowa- 

nie tych wzorów, biorąc do obliczeń większą liczbę pochodnych z sze- 

regu Taylora przy przechodzeniu od równania (23) do równania (26). 


Wykazał on poprzez szereg obliczeń numerycznych [33], że liczba półek 
teoretycznych obliczona z: 


(26) 


(26a) 


(23) 


. 
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(т) 
Ni (У) 


wymaga stosowania poprawek dla większych obciążeń kolumny. I tak: 
(a) dla р<0,75 


= NT Pp Б 2 211 3 
Нал [2+ 3 Г45Р 1 9835 Р (27) 
(b) dla p < 1,6: 
ШУР Eh T a > 
Pra [+ 3 Т 45Р t 1—0,565p 


(c) dla obciążeń większych, tj. dla р = 1,6, całkowitą liczbę półek można 
obliczyć z: 


(28) 


478, 
(М,—У,)? 


(gdzie znaczenie poszczególnych zmiennych jest pokazane па гус. 7). 


ВЕ (29) 


krzywa 
- przebicia 


V2 Vo,s V4 


Ryc. 7. Graficzne przedstawienie zmiennych równania (29) 


Teoria Mayera-Tompkinsa [2, 3]. Założono, że rozwa- 
żana kolumna podzielona jest na szereg hipotetycznych odcinków, z któ- 
rych każdy odpowiada jednej półce teoretycznej i zawiera 5% ml roz- 
tworu oraz òm g jonitu. Ponadto roztwór przepływa skokowo, porcjami 
ро ёо ті], przy czym nowy skok zachodzi każdorazowo po ustaleniu się 
równowagi. 

W stanie równowagi na półce n przy porcji r roztworu 5% ilość sub- 
stancji zawartej w roztworze określona jest przez: 


== $®С, 
w jonicie zaś przez: 


© 
Sr,n = ûm di 
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Stąd współczynnik podziału dla całej półki można określić przez 
Ku Srm 


Lyn 


Równanie bilansu dla takiego modelu można napisać w formie [2]: 


Lr,n + Srn = Lr, п—1 a Sr—in 
Następnie należy dokonać takich przekształceń, aby otrzymać równanie 
krzywej elucji w formie: 

| Lyn = 9 (п,т,К м) 
W tym celu przyjmiemy, że brana pod uwagę substancja zaobsorbowana 
jest na wierzchołku kolumny, gdzie tworzy wąskie pasmo, które umow- 
nie określono jako półka 0, zaś ilość substancji na niej zawartej jako 1. 
W takiej sytuacji przed rozpoczęciem rozwijania całkowita ilość zaabsor- 
bowanej substancji znajduje się w fazie jonitu na półce 0, co oznacza, 
że 5: ,=11 І, = 0. Równanie bilansu będzie miało wtedy postać: 


L, b 8,,F0+1 
korzystając z definicji współczynnika podziału mamy: 


8,0 
U = к 
І, ża 
stąd: 
М Ki 1 
w TEREN IFE y 
podobnie: 


1, +8. И. 9 oraz : Szo =K 
2,0 2,0 1 — Ka 5 Li M 


Śledząc dalej zmiany Г i 5 na pierwszej półce zależnie od ilości” 
przepływającego przez nią roztworu, a następnie zmiany, jakie wywo- 
łuje przechodzenie 5% ml roztworu przez kolumnę, można dojść do wzo- 
rów ogólnych [2]: 

i (+n—1)! К 
мк penna (ү ҮЗ" 
(r—1)!n! (1-Кы) (30) 
и (r+n— 1)! еы 
бэт: (r-l)ln! (1+K,)*” 
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Dla dużych wartości r i n można zastosować przybliżenie Stirlinga, co 
pozwoli pozbyć się kłopotliwych w obliczaniu silni (п! = е "n"V/2zn ) 
Funkcja (30) ma kształt krzywej Gaussa z maksimum w: 


Lmax = Ln,n [2:7 n: Km (1 + Ku)] + 


Stąd po pewnych przekształceniach otrzymamy: 


Км т? А 
не (е ЧИ, 
lub dla dużych Км: 
czę KFS 
к=з | (31b) 


Teoria Glueckaufa dla rugowania [5]. Załóżmy istnie- 
nie urojonego ruchu jonitu z taką prędkością — w, aby można było 
przyjąć, iż izochrona pozostaje w jednym miejscu na kolumnie. Oznacza 
to, że ruch jonitu jest równy co do bezwzględnej wartości szybkości 
migracji stacjonarnego pasma zaabsorbowanej substancji, lecz odbywa się 
w przeciwnym kierunku. Wtedy przez przekrój kolumny (ryc. 8) w punk- 
cie 1 opuści warstewkę dl jonit o stężeniu 9; (1) i równocześnie przez 
powierzchnię w punkcie 1+1 wpłynie do niej jonit o stężeniu q; (l+ dl). 


l+dl 


Ryc. 8. Schemat kolumny ze stacjonarną izochroną; A — kierunek urojonego prze- 
pływu jonitu z prędkością -w, B — kierunek przepływu roztworu 
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W jednostce czasu, wskutek przepływu jonitu, zawartość substancji 


w warstewce dl ulegnie zmianie o: 1 
0q,| _ 
-.| 3t). SW [9:05 dl) — q,(0)] 


W tym czasie taka sama ilość substancji musi ulec reakcji wymiany, 
ponieważ pasmo, zgodnie z przyjętymi założeniami, ma pozostać w tym · 
samym miejscu na kolumnie. Stąd: 
5a) 
"a=; 32 
к, 5 


$ (2) т. 
dl 


Jak już wzmiankowano, szybkość wymiany jonowej jest zdetermino- 
wana przez dyfuzję. Z tego powodu dla procesów wymiany, zdetermino- 
wnaych przez dyfuzję w warstewce adhezyjnej, można prawą stronę 
równania (32) zastąpić przez równanie (9c). Następnie, korzystając z de- 
finicji współczynnika separacji (К), należy znaleźć relację typu 
qi= 9 (Ks, С;), która pozwoli z równania (32), po uwzględnieniu w nim 
dyfuzji, wyrugować zmienną q;. Po scałkowaniu tego równania oraz nor- 
malizacji przeprowadzonej tak, aby początek nowego układu współrzęd- 
nych M się zawsze w środku izochrony (ryc. 9), otrzymamy: 


_qlC_ _ 0,266 тоЕ i K 1 1 1 
U (== [== In — а], 
ih D,(1+70r.F | | Ks—1 PoR Ф Mdr 


-15 -10 -5 о 5 10 U/zi/ 


Ryc. 9. Kształt izochrony obliczonej z równania (33) dla ЕНТРр + EHTPo = 4 ст, 
ЕНТР; + ЕТТР, = 2 cm, Ks = 10. Współrzędna Г(2) została znormalizowana w taki 
sposób, że środek ciężkości izochrony ma współrzędną (2) = 0. [6] 


Glueckauf w miejsce wyrażenia wziętego w nawias kwadratowy wstawił 
wartość Ну (patrz równanie 19). Zważywszy jednak, że wyrażenia te 
różnią się o: 
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qo/C 
qo/C+f 


a to przy f —),4 jest prawie równe jedności. Dlatego dla zachowania 
jednolitości teorii wygodniej jest przyjąć wartości Glueckaufa [5] 
bądź Helffericha [6]. Ostatni uważa, że uzyska się lepszą adekwat- 
ność teorii, jeśli do Н; doda się Но, co pozwala na uwzględnienie wpływu 
rozmiarów ziarna jonitu na kształt krzywej przebicia. 

Podobnych operacji obliczeniowych dokonuje się dla procesów zde- 
terminowanych przez dyfuzję wewnątrzcząsteczkową, z tą różnicą, że tym 
razem prawą stronę równania (32) zastępujemy przez równanie (9b). 
W wyniku rozwiązania (po uwzględnieniu poprawek Glueckaufa i Helffe- 
richa) otrzymamy: 


Hp+Ho (1 орет. РЕИС SB 
e=- zg а TDM Г, 


lub ogólnie: 
НҮНӨ: ryad "Эна <<< 1 
(2) = ә (Ea тігі) + 


ps 


2 КУС ua Ks-1 Miz) (33) 

Teoria Speddinga-Tremillona [8, 16, 34—36]. Rozpatrz- 
my wyodrębnioną na kolumnie warstewkę, w której leży izochrona, przy 
czym załóżmy, że rozwijanie trwało na tyle długo, iż zdołał wykształcić 
się stan stacjonarny oraz, że przepływ jonitu odbywa się w kierunku 
przeciwnym do przepływu roztworu. Wtedy równanie bilansu napisać 
można w postaci [8, 16]: 


BV Си + бт ати = V Cin) + дт Ят) 


Oznaczmy sąsiadujące ze sobą jony odpowiednio: przez 1 — jony 
silniej absorbujące się i przez 2 — jony słabiej absorbujące się. Wtedy 
na krańcach izochrony będziemy mieli (ryc. 8) powyżej izochrony: 91= qo; 
С, = Co, а poniżej izochrony qı = С, 0, 

Po uwzględnieniu ostatnich relacji w równaniu bilansu oraz dodatko- 
wych przekształceniach [8, 16] można otrzymać zależność: 


o A 
Чиа Co” Cim 


Dla całej kolumny można napisać takich równań n, tj. tyle, ile półek liczy 
dana kolumna. Jeżeli z każdego równania wyrugujemy qı w oparciu 


82 Włodzimierz Hubicki, Lucjan Pawłowski 


o definicję współczynnika separacji i następnie pomnożymy je przez 
siebie stronami, otrzymamy: 


ulet Cim _ Си _ 
(7 [et а-а w 


Jeżeli przyjmiemy założenie, że na półce p leży środek ciężkości izo- 
chrony, tj. punkt o stężeniu С, = 0,5 Co, otrzymamy [8, 16]: 


Cin EN 
Е S Ea 4 
(С. С.) (= | mo 


gdzie п w tym przypadku oznacza liczbę półek, o jaką .oddalony jest 
od środka ciężkości izochrony punkt o stężeniu C,=0,5 Со. Po zlogaryt- 
mowaniu równania (34) otrzymamy: 


$ 1 
In, - INR, =1n Ks: (ат) \ (35) 


W podobny sposób można postąpić z równaniem (34а). 
Wykreślając (In К, —Шш В.) przy l otrzymamy linię prostą, której kąt 
nachylenia do osi 1 będzie określony przez maš co pozwala na ustale- 


eksperymentalnej wartości EHTP w prosty sposób. 


OKREŚLANIE DŁUGOŚCI KOLUMNY DLA PRZEPROWADZENIA ROZDZIAŁU 
Z UZYSKANIEM PRODUKTÓW O ŻĄDANEJ CZYSTOŚCI [5, 6, 29, 33] 


Rozważmy dwa sąsiednie pasma z zaabsorbowaną substancją 1 i 2. 
(ryc. 10). Zanieczyszczenia bezwzględne określone są przez 


Am, al Ат; 


ыл эдин тщ (substancja 2 przez 1) 


теди Ат, Ат, 
1 = т, — Ama z Mı 


(substancja 1 przez 2) 


Dla 1: = т po uwzględnieniu równania (23) i założeniu, że ЕНТР,= ЕНТР,, 
otrzymamy: 
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С,=0,368 с, тах 


Ryc. 10. Schemat krzywej elucji 


Ostatnią zależność ilustruje ryc. 11. 


N wupu 


2 
о > 
+ u 


ub шога м 


z 

——> 
ё > 
w 


= гал 4 


Ryc. 11. Diagram dla małych obciążeń kolumny [6] 
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W rzeczywistości ЕНТР, 52 ЕНТР. i wtedy [29]: 


m; + mą _ + AZ 
"BRT, r. asy (У. +7) 
gdzie N jest średnią arytmetyczną z liczby półek teoretycznych dla - 
obydwu substancji. W tym przypadku wielkością odczytaną z diagramu 
(ryc. 11) jest N=(N;+N;):2, stąd optymalna długość kolumny wyniesie 
I=N *EHTP.. Stosowanie diagramu z гус. 11 jest ograniczone przez 
warunek (26a), ponieważ obliczenia, do niego prowadzące, zostały oparte 
na tych samych przybliżeniach. 

Dla większych obciążeń kolumny a, obliczanie optymalnej długości 
kolumny jest bardziej złożone. Przeprowadził je Glueckauf [38]. 
Dla przykładu podał on diagram (ryc. 12), w którym określił maksymalne 


| 
Гы 
ca 
a 
ка 
Мы „| 

е9 Б-ы 
РУ N 


РН 
5 Ех ЫМУ НРУ 


А Е 
чу чи 


У 


DY 
ГУ 
RRT RE 


п 
і 
[М 
a 
DN 
LIP 


1073 2. 34 68102 2 34 SZA 34 
ле lub © 
Ryc. 12. Diagram dla dużych obciążeń kolumny [33] 
obciążenie kolumny №, pozwalające przy danej długości kolumny М 
i współczynniku separacji К, otrzymać produkty o czystości 1=99,9%. 


UWAGI KOŃCOWE 


Z przytoczonych danych wynika, że kształt krzywej przebicia zależy 
od efektywnej wysokości teoretycznej półki (EHTP) i od współczynnika 


Zagadnienie optymalizacji jonowymiennego rozdziału pierwiastków... 85 


separacji (К;). Jednak wysokość półki może być również określona przez 
UN, gdzie N w zależności od konkretnej sytuacji znajdujemy z równań 
(25—29, 31a, 31b, 34, 34a), Należy podkreślić, że wartość EHTP otrzy- 
mana na drodze teoretycznej (równania 13—16, 18—21) zwykle dość 
znacznie różni się od wartości empirycznej. Mimo tych różnic są one 
jednak bardzo użyteczne przy badaniu wpływu takich parametrów na 
jakość rozdziału, jak: współczynnik separacji, szybkość przepływu 
eluentu, rozmiary ziarna jonitu i temperatura. 

Wysokość półki teoretycznej, określona na drodze eksperymentalnej, 
jest również miarą sprawności kolumny i z tego względu była wielo- 
krotnie brana pod uwagę [2, 5, 16, 17, 19, 22, 23, 25, 26, 27, 28, 34—44]. 


ZAGADNIENIE OPTYMALIZACJI ROZDZIAŁU Eu OD Sm 


W jonowymiennych rozdziałach pierwiastków ziem rzadkich stosuje 
się zasadniczo dwie techniki: elucyjną i rugowania. W wielu przypad- 
kach spotyka się jednak w piśmiennictwie z pewnym pomieszaniem Sensu 
tych pojęć. Według Helffericha [6] rugowanie następuje wtedy, 
gdy jony roztworu rozwijającego mają większe powinowactwo do jonitu 
od jonów zawartych w paśmie rozwijanym, natomiast elucja — wówczas 
gdy jony roztworu rozwijającego mają mniejsze powinowactwo do jonitu 
od jonów zawartych w paśmie rozwijanym. Zasadnicza różnica pomiędzy 
tymi dwoma technikami polega na tym, że w przypadku rugowania 
pasma rozdzielanych substancji ściśle przylegają do siebie w czasie ich 
wędrówki wzdłuż kolumny, natomiast w przypadku elucji pomiędzy nie 
wchodzą jony z roztworu rozwijającego, co daje możność całkowitego 
rozsunięcia pasm; stosowanie elucji w celach preparatywnych jest jednak 
niecelowe, gdyż otrzymany tą techniką wyciek posiada małe stężenie. 

Do celów preparatywnych (przy pracy z dużymi ilościami), gdzie 
wydzielenie całkowitej ilości określonych jonów nie jest konieczne, lepiej 
nadaje się technika rugowania. Jakkolwiek całkowite rozsunięcie pasm 
nie jest możliwe, to jednak przy dużych obciążeniach kolumny część, 
w której pasma na siebie zachodzą, stanowi tylko niewielki procent 
całości. Tak więc zawracanie frakcji z pogranicza pasm do ponownego 
rozdziału pozostaje bez większego uszczerbku dla wydajności. 

Elucja. Omówione teorie Mayera i Tompkinsa oraz wyszczegól- 
nione wyżej niektóre teorie Glueckaufa opracowane zostały pod kątem 
zastosowania ich przy interpretacji elucyjnych rozdziałów jonowymien- 
nych. Celem zilustrowania relacji odnoszących się do elucji przeanalizo- 
wano rozdział europu od samaru w temperaturze pokojowej па Amberli- 
cie IR-120 o rozmiarach ziarna jonitu т,= 5 · 10-3 ст. Współczynniki dy- 
fuzji przyjęto z pracy Cornisha [26]. Stosując równania (13—16), 
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badano wpływ KG na wartość EHTP. Okazało się, że wartość EHTP 
maleje wraz ze wzrostem Kg, zgodnie z ogólnymi zasadami, przy czym 
poczynając od Kg=16 obserwuje się spłaszczenie krzywej i dalej już 
w małym stopniu wartość EHTP zależy od Kg. Z tego wynika, że opty- 
malną wartością współczynnika podziału będzie 16. 

Stosując jako eluent EDTA można tak dobrać jego stężenie i pH, 
że będzie możliwe osiągnięcie Kg o założonej z góry wartości, Współczyn- 
nik separacji dla pary Eu-Sm ma wartość 1,4 [35]. Biorąc Кс=13,3 
i Kg=18,6 spełni się wyżej wymienione warunki (К,, Ко). 

Przy założeniu, że mguFmMsm w oparciu o diagram z ryc. 11 znale- 
ziono odpowiednie liczby półek dla uzyskania rozdziału z otrzymaniem 
produktów o czystości 99,0 i 99,99% (patrz tab. 1). 


Tab. 1. Teoretyczne zależności dla jonowymiennego rozdziału Eu od Sm przy użyciu 
EDTA jako eluentu (pole przekroju rozważanej kolumny 10 cm?) 


e = Czysty ЕНТРЕ, +EHTPSm — C на. 
Е алын > 1 VEu wymycia 
produkt| N 2 
cm3 зек.—1 cm—2 К ст | cm? piku 


(godz.) 


99,0. | 180 0,256 46,1 | 6316 35,1 
0,005 99,9 380 0,256 97,3 | 13330 74,1 
99,99 | 510 0,256 130,6 | 17892 99,4 
99,0 180 0,504 ` 90,7 | 12426 34,5 
0,01 99,9 380 0,504 191,5 | 26236 72,9 
99,99 | 510 0,504 · | 257,0 | 35209 97,8 
99,0 180 1,000 180 | 24660 34,25 
0,02 99,9 380 1,000 380 | 52060 72,3 
99,99 | 510 1,000 510 | 69870 97,0 
99,0 180 1,992 358,6 |49128 34,1 
"0,04 99,9 380 1,992 757 [103709 72,0 
99,99 | 510 1,992 1015 |139178 96,6 


Interesująco przedstawia się zależność EHTP od F. Okazuje się, że 
szybkość przepływu eluentu w znaczny sposób wpływa na wartość EHTP. 
Jak wynika z danych przedstawionych w tab. 1, najlepsze są, jak się 
wydaje, małe szybkości przepływu. Wzrost F przyczynia się do wzrostu 
wartości EHTP, a tym samym do konieczności stosowania coraz dłuższych 
kolumn. To z kolei powoduje wzrost zużycia eluentu i jonitu przy mini- 
malnym skróceniu czasu rozdziału (patrz tab. 1). 

Przeanalizowane zależności nie są ilościowe, ponieważ wartość EHTP, 
obliczona na drodze teoretycznej, nie odpowiada w pełni wartości ekspe- 
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rymentalnej. Głównym celem przeprowadzonej w tej pracy dyskusji 
było zilustrowanie samej metody. Przy doborze optymalnej długości ko- 
lumny należy raczej korzystać z wartości EHTP określonej na drodze 
eksperymentalnej. 

Z przeprowadzonych obliczeń wynika, że otrzymane w oparciu o rów- 
nanie (31а) dane (teoria Mayera i Tompkinsa) różnią się od uzyskanych 
z równań (25) i (26) (teoria Glueckaufa), nawet o kilkadziesiąt procent 
przy małych współczynnikach podziału, np. dla pary Eu-Sm przy współ- 
czynniku podziału równym 4 różnica wynosi ok, 20%. Wydaje się, że 
w dyskusji jonowymiennych rozdziałów pierwiastków ziem rzadkich 
stosowanie teorii Mayera i Tompkinsa mija się z celem z uwagi na to, 
że są do dyspozycji równania Glueckaufa, bardziej adekwatne do danych 
eksperymentalnych. Należy jednak zwrócić uwagę na inną zaletę teorii 
Mayera i Tompkinsa — tłumaczy ona mianowicie bardzo poglądowo 
istotę półki teoretycznej 

Rugowanie. Cechą charakterystyczną rugowania jest to, że jony 
rozwijające mają większe powinowactwo do jonitu od jonów zawartych 
w paśmie rozwijanym, co uniemożliwia im wnikanie w głąb pasma jonów 
rozdzielanych. У 

W przypadku kompleksów z EDTA można poprzez zmianę pH i stę- 
żenia tak daleko zredukować powinowactwo pierwiastków ziem rzadkich 
do jonitu, że w konsekwencji od elucji przejdzie się do rugowania. 
Ze względu na stechiometryczny charakter wymiany jonowej stężenie 
pierwiastków ziem rzadkich w roztworze na kolumnie jest równe stężeniu 
roztworu rozwijającego. Te zaś zawsze można zmieniać zgodnie z po- 
trzebą. 

Przeprowadzona dyskusja równań (18—21) dla pary Eu-Sm wyka- 
zała, że EHTP przy rugowaniu zależy w taki sam sposób od F, то, Di, Dri 
jak w procesie elucji; oraz dodatkowo od stężenia roztworu rozwijają- 
cego, gdzie maleje ze zmniejszaniem się jego wartości — początkowo 
dość gwałtownie, a w miarę wzrostu q,/C coraz wolniej. Obniżanie war- 
tości qo/C poniżej 16 jest, z teoretycznego punktu widzenia, nieopłacalne. 
Z drugiej strony należy liczyć się z faktem, że zmniejszanie stężenia 
roztworu pozostaje w prostej zależności ze stężeniem produktu w wy- 
cieku. Stąd jego nadmierne rozcieńczenie zmusza do obróbki dużych 
objętości wycieku, co jest zjawiskiem niekorzystnym. Posługując się 
określonymi przez równania (18—21) wartościami EHTP,, EHTP, i EHTP; 
można z równania (33) wykreślić teoretyczną izochronę. 

Wykreśloną teoretyczną izochronę dla EHTP;= 0,372 cm, EHTP,= 
= 0,372 cm, ЕНТР,= 0,124 cm i ЕНТР.=0,008 cm przebadano przy użyciu 
równania (35). Okazało się, że wyniki pozostają w dobrej ze sobą zgod- 
ności, gdyż obliczona przy pomocy równania (35) (teoria Speddinga- 
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-Trómillona) wartość ЕНТР=0,51 cm jest prawie równa wartości, jaką 
wstawiono do równania (33) (teoria Glueckaufa), czyli ЕНТР= [ЕНТР,+ 
+ EHTP,+ EHTP;] = 0,504 cm. 


* * 


Przeprowadzone w niniejszej pracy rozważania teoretyczne są pod- 
stawą do interpretacji danych eksperymentalnych, uzyskanych w rozdzia- 
łach mieszanin pierwiastków ziem rzadkich metodą jonowymienną przy 
zastosowaniu różnych pod względem chemicznym eluentów. 

Konfrontacja danych teoretycznych i eksperymentalnych będzie ce- 
lem prac następnych. 


PIŚMIENNICTWO 


Martin i Synge: Biochem. J., 35, 1385 (1941). 

Mayer S. W, Tompkins Е. R.: J. Am. Chem. Soc., 69, 2866 (1947). 

Tompkins Е. В.: Г. Chem. Education, 26, 32, 92 (1949). 

Glueckauf E.: Trans. Faraday Soc., 51, 34 (1955). 

. Glueckauf E.: Ion Exchange and its Applications. Society of Chemical In- 
dustry, London 1955, s. 34. 

6. Helfferich F.: Ion Exchange, McHill, New York 1962. 

7. Іпс2ёау J.: Analytical Applications of Топ Exchangers. Pergamon Press, Lon- 
don 1966. 

8. Trómillon B.: Sóparations par les résines ćchangeuses dions. Paris 1965 
(wydanie rosyjskie: Moskwa 1967; wydanie polskie: Jonity w procesach roz- 
dzielczych. PWN, Warszawa 1970). 

9. Boyd R. Е, Adamson A. M, Mayers L. R. Jr.: J. Am. Chem. Soc., 69, 
2836 (1947). 

10. Bird R.B.: Theory of Diffusion. New York 1959. 

11. Thomas H.C.: J. Am. Chem. Soc., 66, 1664 (1944). 

12. Vermeulen T.: Ind. Eng. Chem., 45, 1669 (1953). 

13. Vermeulen T., Hiester М. К.: J. Chem. Phys., 22, 96 (1954). 

14. Oppler A., Hiester М. K.: Tables for Predicting the Performance of Fixed 
Bed Ion Exchange. Stanford Research Institute 1954. 

15. Glueckauf E.: Trans. Faraday Soc., 51, 1540 (1955). 

16. Trómillon B.: Bull. Soc. Chim. France, 508 (1958). 

17. Trémillon B.: Bull. Soc. Chim. France, 281 (1958). 

18. Jouy A. D,Coursier J.: Bull. Soc. Chim. France, 323 (1958). 

19. Cousier J, Hurre J.: Anal. Chim. Acta, 18, 272 (1958). 

20. Kambara T.: Japan Analyst, 19, 234 (1969). 

21. Cornet C, Cousier J, Huerre J.: Anal. Chim. Acta, 19, 259 (1958). 

22. Tróćmillon B.: Bull. Soc. Chim. France, 1621 (1958). 

23. Trómillon B.: Ann. Chim., 4, 1055 (1959). 

24. Raczyński W. W.: Zarys teorii dynamiki sorpcji i бырай PWNT, 
Warszawa 1966. 

25. Cornish F. W.: Analyst, 83, 634 (1958). 

26. Cornish Е. W.: Canad. J. Chem., 34, 1471 (1556). 


т до № н 


Проблема оптимализации ионообменного разделения... 89 


27. Dybczyński В.: J. Chromatogr., 31, 155 (1967). 

28. Dybczyński R.: J. Chromatogr., 14, 79 (1964). 

29. Glueckauf E.: Trans. Faraday Soc., 51, 34 (1955). 

30. Freund J. E.: Podstawy nowoczesnej statystyki, PWE, Warszawa 1968. 

31. Bronsztejn J, Siemiendiejew K.: Matematyka — poradnik ency- 
klopedyczny, PWN, Warszawa 1968, str. 94. 

32. Handbook of Chemistry and Physics, str. 204, Chemical Rubber Publishing Co. 

33. Glueckauf E.: Trans. Faraday Soc., 60, 729 (1964). 

34. Tremillon B.: Bull. Soc. Chim. France, 502 (1968). 

35. Powell Ј. Е, Spedding Е. Н, Svec Н. J.: Chem. Eng. Prog. Sym- 
posium Series 55, 101 (1959). 

36. Spedding Е. Н, Powell J. E, Svec H. J.: J. Am. Chem. Soc. 77, 6125 
(1955). 

37. Dybczyński R.: Raport 1115 /УПИС IBJ, Warszawa 1969. 

38. Kambara T.: Japan Analyst, 21, 1536 (1963). 

*39. Marciny E.: Duyckaerts, J. Chromatogr. 22, 164 (1966). 

40. Mann С. K.: Anal. Chem., 32, 67 (1960). 

41. Mann С. K, Swanson С. L.: Anal. Chem., 33, 459 (1961). 

42. Powell J. E, Burkholder H. R.:J. Chromatogr., 29, 210 (1967). 

43. Inczódy J.: Ion Exchange in the Process Industries. Society of Chemical 
Industry, London 1969. 


44, Powell J. E, Spedding F. H.: Trans. Metal. Soc. AIME, 215, 457 (1959). 


PE3IOME 


В работе дается характеристика влияния основных параметров 
(т.е. коэффициентов внутримолекулярной диффузии и в адгезион- 
ной пленке, размеров зерен ионита, скорости протока раствора и ко- 
эффициентов распределения) на величину эквивалентной высоты 
теоретической тарелки, которая может служить мерой разделитель- 
ной способности ионообменной колонки. 

Хотя не было найдено согласия между теоретически расчитанной 
величиной ВЭТТ и экспериментально определенной величиной, тем 
не менее первую из них можно считать хорошей качественной ре- 
презентацией влияния вышеуказанных параметров на разделитель- 
ную способность колонки и поэтому она может применяться для под- 
бора их оптимальных величин. 

В работе описываются методы Мейера и Томпкинса, Глюкауфа 
и Спеддинга-Тремийона определения числа теоретических тарелок на 
основе экспериментальных данных. Эти методы позволяют косвенно 
определить эмпирическую величину ВЭТТ, которая является полез- 
ной для точной характеристики разделительной эффективности ко- 
лонки, а совместно с коэффициентом разделения может служить 
для определения оптимальной длины колонки. для заданного разде- 
ления. 
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твором ЭДТА методом элюции показали, ч?о теория Мейера-Томпкин- 
са дает результаты, отличающиеся на несколько десятков процентов 
от данных, расчитанных по уравнениям Глюкауфа. Эти разницы 
уменьшаются по мере роста коэффициентов распределения. 

Для разделения той же смеси этой пары элементов вытеснитель- 
ным методом величина ВЭТТ, расчитанная по теории Спеддинга-Тре- 
мийона, почти идентична с величиной полученной по уравнениям 
Глюкауфа. 


SUMMARY 


The paper presents the effect of fundamental parameters (i.e. coeffi- 
cients of particle diffusion and film diffusion, particle size, flow rate 
of solution and distribution coefficients) on the equivalent height of 
a theoretical plate (ЕНТР) which can be considered a measure of the 
effectiveness of ion exchange columns. { 

Although the theoretically calculated and experimental values of 
EHTP are not in an exact agreement, the former can be considered 
a good qualitative representation of the effect of the above mentioned 
parameters on the separation ability of the column and thus it can be 
used for the selection of their optimum values. 

The methods elaborated by Mayer and Tompkins, Glueckauf and 
Spedding-Trćmillon for the determination of a number of theoretical plates 
on the basis of experimental data have been discussed. The methods 
allow to determine indirectly the empirical value of EHTP which can be 
useful to the precise characterization of the effectiveness of a column, 
and which in combination with the separation factor can be employed 
for the determination of optimum column length of a given separation. 

Theoretical calculations of the separation of the pair Sm-Eu by elution 
with EDTA have shown that the results of Mayer and Tompkins’ theory 
differ in several ten per cent from data following from Glueckauf's 
equations. The differences decrease with the increasing values of distribu- 
tion coefficients. 

The value of EHTP calculated from Spedding-Trómillon's theory 
for the separation of the above mentioned pair is almost identical with 
that obtained from Glueckauf's equations. 
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